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OO ' Abstract 
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. A relativistic wave equation for bound states of two fermions with arbitrary masses 

which are exposed to a magnetic field is derived from quantum electrodynamics. The inter- 
action kernels are based upon the generalized invariant M. -matrices for inter-fermion and 
fermion-field interactions. As an application we calculate the energy corrections in a weak 
homogeneous B field to obtain the Zeeman splitting of the hyperfine structure (HFS) and 
. g-factors in the lowest order (i.e. to 0(a 4 )). Lande ^-factors are presented for several of 

the first excited states of hydrogen, muonium, and muonic-hydrogen. 

^3 ' 
p 

1. Introduction 

Ph. 

The relativistic treatment of energy levels of two-fermion atomic systems (including 
j> [ atomic hydrogen, hydrogen-like ions, helium-3 ion, muonium, muonic-hydrogen), as well 

as their fine structure (FS) and hyperfine structure (HFS) in an external uniform magnetic 
field (Zeeman effect), is an important problem. The theoretical knowledge of energy spectra 
and transition frequencies provides a test of two-body bound-state QED [1]. One can then 
obtain information about the character of the coupling in the system, the gyromagnetic ra- 
tios of the bound particles, the magnetic moments [2], the mass ratio [2-5], and fundamental 
physical constants such as the Rydberg constant Roo, and the fine structure constant a [6]. 
The Zeeman effect in the HFS can be used as a diagnostic tool for solar photospheric mag- 
netic fields [7], fusion research and plasma physics, where the magnetic field is applied to 
control the shape and position of the plasma [8]. 

In the lowest-order approximation the linearly dependent part of the energy splitting for 
a two-fermion system placed in a weak static magnetic field B can be written as [1,9-11] 

^ E F,lnj,hAsuI = iVBigi + /iB2#2) Bm F , (1) 

where F, mj, ji, £, s\, I are quantum numbers, which characterize the system: si and / are 
the spins of the first and second particle respectively, I and j\ represent the orbital and total 
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angular momentum quantum numbers of the first particle. The total angular momentum 
of the system is denoted by the quantum number F — j\ + I, ji + I — 1, \ji — I\. The 
projection of the total angular momentum on the B direction is mj? = —F, — F+l, ...F—l, F. 
The "Bohr magnetons" for the two particles are defined as /j,bi = Qih/2m 1 c, and /j,b2 — 
—Q2h/2m 2 c, where Qi, Q 2 > 0). Usually, in our notation m 1 and m 2 correspond to the light 
and heavy particle respectively. Assuming that the energy-level splitting (1) is much smaller 
then the HFS splitting, AE ext « AE HFS , the Lande (g-) factors g\ and g 2 take the form 
[9-11] 

= F(F + l) + j 1 (,- 1 + l)-/(J + l) 
yi yji 2F(F + 1) v ; 

where 

0* - i + - 1) iiC7i + i) + M3 1 + i)-^(^ + i) 

^ 1 ^ 2j! (j! + 1) ' 1 ^ 

and 

„ „ F(F+l)- 3l ( Jl + l) + /(/ + l) 

^ 2 -fe 2F(F + 1) • (4) 

Here y" si and y- S2 are the intrinsic spin magnetic moments of the constituent particles. Ac- 
cording to the Dirac theory a free particle at rest has g s = 2. In QED g s is corrected by the 
anomaly, which to lowest order is given by the Schwinger correction. For bound particles 
the intrinsic moment can be expressed as 

*m = 2 + + A ^> (5) 

where the terms Ag REL , Ag QED represent the relativistic [9,12,13], and QED corrections 
respectively (cf. the review [14]). There is also an additional higher-order contribution to 
(1), Ag±2 fJ'Bi^BrriF, which is caused by the hyperfine structure (HFS) [15]. 

The g-factors (2) and (4) are not symmetrical, because they were obtained under the 
assumption that the orbital motion of the heavy particle can be neglected. In hydrogen the 
nucleus contributes a fraction of m\j (mi + m 2 ) ~ 5 x 1CT 4 to the orbital angular momentum, 
while for muonic hydrogen this fraction is ~ 0.1. The relativistic and QED corrections in 
(5) can be comparable with the orbital angular momentum effects of the heavy particle. 
Recent high-precision measurements of the ^-factor in hydrogen-like systems have reached 
an accuracy of about 5 x 1CT 9 [16,17]. Thus, it is desirable to obtain a more general result 
for the g- factor in order to overcome the shortcomings of Eqs. (2,4). It will be shown that 
this is particularly important for excited states. 

In this work we present an analysis of the HFS of a two-fermion system in an external 
magnetic field based upon a reformulation of QED and the variational Hamiltonian formalism 
developed earlier [18-20]. A relativistic two-fermion wave equation for arbitrary fermion 
masses is, thus, derived from first principles. A solution of this equation permits, in principle, 
to obtain all QED energy corrections to any order of the coupling constant [18]. In the present 
paper we extend the method to derive the integral wave equation in momentum space for 
the case where a uniform weak magnetic field is present. We calculate the Zeeman splitting 
of the HFS energy levels to O (a 4 ) for all quantum states and unrestricted values for the 
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fermion masses. We obtain a novel result for the g-factor, Eqs. (38-41), and demonstrate 
that it coincides with Eqs. (2-4) in the case of m 2 >> mi, as long as the intrinsic moment 
of mi is restricted to the Dirac value g s = 2. 

The modification of the wave equations due to the external magnetic field is presented 
in Section 2. In Section 3 we provide the classification of the quantum states, and a partial- 
wave decomposition of the momentum-space equations. Section 4 contains expressions for 
the Zeeman energy splittings of the HFS levels, and the (/-factor results. Numerical values 
for the Lande factors are compared with data from Eqs. (2,4) for various excited states 
of hydrogen, muonium and muonic hydrogen. In most expressions we use natural units 
h = c=l. 



2. Bound-State Variational Wave Equation 

For two-fermion systems without external fields wave equations were derived in [18-19] on 
the basis of a modified QED Lagrangian [21-22]. With this Lagrangian a simple Fock-space 
trial state 

Atrial) =Y,j d' Vl d' V2 F slS2 { Vl ^ 2 )bl iS Dl 2S2 |0> , (6) 

sufficed to obtain HFS levels correct to fourth order in the fine-structure constant. Here 
6 qiSl and -D q2S2 are creation operators for a free fermion of mass m\ and an (anti-)fermion 
of mass m 2 respectively, and |0) is the trial vacuum state such that 6 qiSl |0) = D^ 2S2 |0) = 0. 

As discussed in section 3 below, the four adjustable functions F SlS2 must be chosen so 
that the trial state (17) is an eigenstate of the relativistic total angular momentum operator, 
its projection, and parity (as well as charge conjugation for the case rrt\ = m 2 such as 
positronium). 

A variational principle is invoked to obtain a momentum-space wave equation for the 
amplitudes [18]: 

= S/ ^ 3 p^ 3 P2K 1 + ^ 2 -^)i 7i slS2 (Pi,P2)^; 1S2 (Pi,p 2 ) (7) 

_ m x m 2 sr^ f d 3 p 1 d 3 p 2 d 3 q 1 d 3 q 2 

( 27r ) 3 <nwi« ' V^WAS 
x ^ia 2 (qi,q 2 ) H) M SlS2Uia2 (pi,p 2 ,qi,q 2 )^* lS2 (Pi,P2), 

where = p\ + m\ and = p\ + m|. The interaction is governed by the generalized 
invariant .M-matrix M. SlS2(71(72 (p 1; p 2 , qi, q 2 ). It has the form 

M { ^ S2>7ia2 (pi, p 2 , qi, q 2 ) = Mf^^ (pi, p 2 , qi, q 2 ) + (Pi' Pa. ^ %) . ( 8 ) 

where M !^ a a (pi, p 2 , qi, q 2 ) is the usual invariant matrix element, corresponding to the 
one-photon exchange Feynman diagram [19-20]. 
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The element A^^* 2Cri(72 represents the interaction with a given external classical field A e * f 
M e s x 1 l 2tTia2 (pi,P2,ci 1 ,q 2 ) (9) 



i (2ix f 



3/2 ( ^^Af( Pl - qi )u ( Pl , Sl ) (-<Q0 t"« ( qi , (Ti) <W 2 



F^pAf (q 2 - p 2 )V (p 2 , <r 2 ) (-<Q 2 ) y V (q 2 , S2 ) 5 



SlO"! 



The Ansatz (6) can not accommodate processes that include the emission or absorption 
of real, physical (as opposed to virtual) photons. Such radiative processes could be included 
by generalizing the trial state. Here we limit ourselves to the form (6), i.e., the effects of 
radiative decay or absorption of radiation are ignored in the present work. 

In order to obtain the Lande factors we evaluate the M.f^ 2(71(T2 matrix (9) in a stationary 
uniform magnetic field B = Bz. The vector potential can be chosen as 

AT (x) = ~yB, AT (x) = l -xB, Ag* (x) = Af* (x) = 0. (10) 

The inverse Fourier transform of the non-zero components yields 

«k) = *!p?,(M^W, <(k) = _w!»^) j(yj(y 

Using the semi-relativistic expansion 

M(pi,Si)7 1 M(qi,ai) = 1^-^P\ X x (Pi - Qi)] + Qi + Pi)i¥><n, ( 12 ) 

M(pi,Si)7 2 M(qi,^i) = 2^^i t^ 1 X ~ qi ^ +qi + Pi)2^i, 

where (p\ = [1 0], y? 2 = [0 l],and (ujp^qA 1 ^ 2 ~ mi, and a similar expansion for anti-particle 
spinors we obtain 

AC' (Pi,p 2 ,qi,q 2 ) (13) 



(2tt) 3/2 / ( Pl - qi )^ (i [?i x ( Pl - qj)] + qi + Pl ) . 

2c ^ + jAf (q 2 - p 2 ) X t 2 (i x (p 2 - q 2 )] + q 2 + p^. Xs2 5 Sl(71 

where x! = [0 1], X2 = — [1 0]> an d J = 1> 2. It is straightforward to show that 

(q^. Af ( Pl - qi ) = - ^^^ L^ ( qi ) 5 ( Pl - Ql ) , (14) 

and 

Af\vi - qi)^ (i [ai x ( Pl - qi )] + qi + Pl ) j <p ai (15) 
= - (2tt) 3/2 B (yl^izip^ + <5 Sl( ji^i« (qi)) ^ 3 (Pi - qi) , 
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where L lz (q x ) is the z-component of the angular momentum operator of the particle with 

mass mi 

L U (Qi) = ~i - ^) ■ (16) 

Taking ip si to be the eigenstates of the spin operator S± z = |<?i 2 , and using a similar 
procedure for the second particle, we obtain 

(Pi,P2,qi,q 2 ) (17) 
(2tt) 3 B ( ^ (MA^i + K^Liz (qi)) 5 S2 aJ 3 (pi - qi) 
~ ^~ \ -£ (2xi 2 S2 Z Xs 2 + S a2S2 L 2z (q 2 )) 5 siai S 3 (p 2 - q 2 ) 

or 

text 



s 1 s 2 <ria 2 (pi,P2,qi,q 2 ) (18) 
( fi B1 (2m ai + L lz (qx)j S 3 (pi - qi) 
-fj,B2 (2ma 2 + L 2z (q 2 )) S 3 (p 2 - q 2 ) 



= -(27r) 3 fi 



where the spin projection quantum numbers rh a can take the values ±1/2. The quantities 
Ubi and fiB2 are the "Bohr magnetons" defined in the previous section. As expected, a 
unit of spin interacts with a magnetic field twice as strongly as a unit of orbital angular 
momentum. 

By going to the next order in the expansion of the invariant M. matrix one can obtain 
self-energy corrections, which lead to divergent loop integrals that have to be cured by 
charge renormalization. The vertex term modifies the Dirac value of the magnetic moment 
by a factor (1 + k), where k is the anomaly (Schwinger correction). This factor can be 
included in our calculation by a replacement 2rh ai and 27710-2 in Eq. (18) by g Sl fh ai and g S2 fn a ^ 
respectively, where g Sia /2 = 1 + fci i2 . The anomaly is the lowest-order QED correction to 
the g factor A^f D = 2k lj2 in Eq. (5). 

3. Partial-wave decomposition and radial wave equations 

The present work is an extension of Ref. [18], in which the partial- wave decomposition of 
the wave equation has been provided. The external magnetic field is treated as a first-order 
perturbation which implies that the quantum labels for the eigenstates do not change. The 
restrictions on the magnetic field strength to justify a perturbative treatment of Eq.(18) are 



B < min 



4 2 4 2 

a m r c a m r c 



H>B1 ^B2 



(19) 



where a = QiQ 2 /^ 1 and m r = m\m 2 j (mi + m 2 ) is the reduced mass. A more explicit 
restriction on B will be presented in Section 4. 
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As outlined in Ref. [18] the trial state (6) is taken to be an eigenstate of total linear 
momentum P, total angular momentum squared J 2 , its projection J 3 , and parity V . It is 
natural to work in the rest frame, where the total linear momentum vanishes. In this frame 
the adjustable functions take the form -F slS2 (pi, P2) = S (pi + P2) F slS2 (p 1 ), where F SlS2 (p 1 ) 
(using p! = p) can be written as 

^ 2 (p) = EE ^ mai32 (p)C 2 (p)- ( 2 °) 

£s 1 s 2 m s 1 s 2 

and Y^ 132 (p) are the usual spherical harmonics. Here and henceforth we will use the 
notation p = |p| etc., while four-vectors will be written as p M . The orbital indices £ slS2 a.nd 
m SlS2 and the radial functions fslll 2 ™ 3132 (p) depend on the spin variables si and s 2 . In the 
rest frame, the operators L\ z (q) and L 2z (q) can be expressed in terms of the orbital angular 
momentum operator, L z (q), of the relative motion: 

L U (q) = — ^— (q) , L 2z (-q) = ^^-L z (q) . (21) 

The substitution of the partial-wave expansion (20) into the rest-frame form of Ansatz 
(6) leads to two categories of relations among the adjustable functions F SlS2 (p): 

(i) The spin-mixed (quasi- singlet and quasi-triplet) states 

In this case we have £ SlS , 2 = I — J, and the general solution under the condition of well- 
defined P, J 2 , J3, and V can be expressed with the help of Dirac T matrices as [18] 

F S1S2 ( P ) = v,rS (p) v -^fj%) + u PS1 rZi (p) v -r»fP(p)- ( 22 ) 

Here fj S9 \p) and ff r \p) are radial functions to be determined. They represent the contri- 
butions of spin-singlet and spin-triplet states, i.e., the total spin is not conserved in general. 

(ii) The t-mixed triplet states 

These states occur for £ S1S2 = t — J =f 1. Their radial decomposition can be written as 

F S1S2 (p) = u pSl T J m -l 2 (p) V^sjj^p) + u psi T J m + s l 2 (p) V_ pS2 f J+1 (p). (23) 

The system in these states is characterized by J, m,j, and V = (— 1) J , and £ is not a good 
quantum number. The two radial functions fj-i(p) and fj+i{p) correspond to the cases 
£ = J — 1 and £ = J + 1. Mixing of this type occurs only for principal quantum number 
n > 3. 

From the variational method we obtain a system of coupled radial equations expressed 
in matrix form as 

( Wp + n p -£)F(p) = ^ f 1 2d 1 K q) F {q) t (24) 

(271) J ^UOpUOqilpilq 
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where cj 2 = p 2 + m 2 and Vll = p 2 + m|, and g = |q| as already mentioned. Here F (p) 
and IK (p, q) are matrices composed of radial functions and kernels respectively. The kernel 
matrix K = K° pe + K ext is made up of one-photon-exchange and external-field parts. Ex- 
plicit expressions for K ope can be found in Ref.[18], while the external-field contributions are 
calculated in this work. 

For the spin-mixed states the two-component Fock-space amplitude is given as 

F(p) = 

The equations imply a mixing of spin and radial variables, and the radial equations are 
usually coupled. We apply a unitary transformation with rotation angle (3 to the spin part 
of function (22) to diagonalize the kernel-matrix. The diagonalization can be carried out for 
arbitrary p and q (cf. Eq. (55) in the Appendix), and defines a new quasi-spin basis 

\si, s 2 , £, s, J, mj) = d \si, s 2 , £,S = 0, J, mj) + C 2 \s 1 ,s 2 ,£, S = 1, J, mj) , (26) 

where £ = J, S is the total spin of the system, and s = for quasi-singlet and s — 1 
for quasi-triplet states. The coefficients used to express the new basis states in terms of 
the previously defined singlet and triplet states are found to be C\ = a/(1 + £) /2, C 2 = 
-a/(1 - /2, for the quasi-singlet states, and d = a/(1 - /2, C 2 = a/(1 +0 /2 for the 
quasi-triplet states. Here the rotation angle (3 has been replaced for convenience according 
to tan 2/? = ^fl^/i. 

The quasi-singlet and quasi-triplet states are both characterized by the same quantum 
numbers J, mj and V = (— 1) J+1 , and they mix the states given in the LS coupling represen- 
tation. The states are labeled for convenience not by the quasi-spin z-projection t 3 = =Fl/2, 
but rather by s = £ 3 + 1/2, which takes on the values of 0,1. In the Appendix the kernels 
for spin-mixed states are given explicitly in order to solve for the angle f3, i.e., to determine 
the ^-values. 

In the limit m 2 » mi the total angular momenta of the first and the second particles 
are j\ — £i± 1/2, j 2 = s 2 = 1/2, where i\ = I. In this case j\ can be used as a good quantum 
number, and the role of the indices s s , s t are played by ji = £± + 1/2 and ji — £\ — 1/2 
respectively. In this case the coefficients C\ and C 2 reduce to C-G coefficients 

C 1>2 = (-l) 1/2+1/2+ ^V(25 + l)(2 Jl + l){ X f I } . (27) 

Note that the one-body limit corresponds to the jij 2 coupling representation, which can not 
be used in the general case of arbitrary masses since ji and j 2 are not independent (they 
are related through the common angular momentum £). For positronium the quasi-states 
become true singlet (C 2 = 0) and triplet (C\ = 0) states with different charge conjugation 
quantum numbers. 



f ( f 9 \p) 
ff'\p) 



(25) 
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We now proceed to calculate the kernels Kf^ n (p, q) associated with the classical external 
field -Aff*. Using Eq. (9) for Mf^ 2aia2 taken in the rest frame, we obtain 



-ext 



g (p, g) - - jT /2y , 2 / ^P' /; q (28) 

iv (mim 2 J 

x Tr 



Q lv /lVV4-*(p - q) ( 7 V + mi ) 7 M (7 A 9a + mi) r» (q) ( 7 A g A - m 2 ) T' m (p) 
-Q 2v ^V 4 ^ t (q " P) (l X qx + mi) r» (q) ( 7 A g A - m 2 ) 7" (7 A $a - m 2 ) T /m (p) 

where g = (a> p , q), and g = — q). The T -matrices correspond to the various J v states. 
The evaluation of these kernels would allow one to obtain all relativistic corrections to the 
g-factor (5), however this is a formidable task. To determine the lowest-order effect it is 
sufficient to use the nonrelativistic limit (q 2 /m 2 « 1). In this case the kernels (28) take 
the form 

^ (p, g) = - {jL jf- / (29) 

/ g^f (p-q)( 7 + /)r( 7 ° + /)r™(q)(7 -/)r /m (p) 
xlr { -Q 2 A^( q - P ) ( 7 ° + /) r» ( q ) ( 7 ° - /) r (7° - 1) r /m (p) 

These are evaluated for a stationary uniform magnetic field (10). The results are given 

separately for the two types of states: 

(i) The spin-mixed states (£ = J, J > 1, V — (— 1) J+1 ) 

In contrast to K*° pe ) (p, q) the kernel matrix K^*) (p, q) is not diagonal in the basis of the 
quasi-singlet \sg q ) and quasi-triplet \tr q ) states, and can be written as 

/CSf } (p, ?) (30) 

/o \3 / Si /Y 1 _ kll ^ m 2 _|_ £fl_ ( \zj _ o |r»i-m 2 | ^ N 
_lf^J_ I mi VV 2J(J+l)y mi+m 2 ^ 2 ^2J(J+1) ^ mi+m 2 ^ y 

2c 1 Sl/Vl l-( \ m-, , ffs 2 / l-£ o |mi-m 2 | ^ 

\ m 2 2J(J+1) j tni+m 2 T 2 ^ 2 J(J+1) ^ m 1 +m 2 ^ 



Bm , 



4f } (p,g) (31) 

In \3 / Qi /Y 1 — l+£ A ™2 I g«i f 1+| I o |rai-m 2 | 

_ if^J_ I 2m lC 2 J(J+1) y mi+m 2 ~T 2 ^2J(J+1) ^ ^ mi+m 2 y ■ 

\ 2m 2 c ^ 2J(J+1) J m 1 +m 2 ~ 2 ^2J(J+1) ~ mi+m 2 



/c<f> (p,(z) = /c 







(27T) 3 


f M 




mi \ 


k ^/J(J+1) 2 


2c 


_Q 2 






y m 2 


\y/j(J+l) 2 



(32) 



2 / mi— m 2 

mi+m 2 



^ (1-^) 



\ 



/ 
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Thus, it couples the system (24). 

(ii) The pure triplet and l-mixed states [I = J =F 1, J > 1, V — (— 1) J ) 

The system (24) can not be decoupled for these states, and the matrix MS ^ (p, q) is not 

diagonal [19]. The magnetic part of the kernel is, however, diagonal 



(p, q) 



(2 



7T 



2c 



roi 


-9i) 


' 1 


" 




m 2 J 





1 



Bm , 



(33) 



All kernels MS ext ^ vanish in the case of equal masses and opposite charges (Q\ = Q 2 ), as 
occurs in the positronium case, where magnetic effects appear only in O (B 2 ) [23]. 



4. HFS to O (a 4 ) order in a magnetic field 

To obtain results for energy levels to O (a 4 ) we solve the radial equations (24) perturba- 
tively using hydrogen-like radial functions (non-relativistic Schrodinger form f^ mj (p)) in 
momentum space [9]. The energy eigenvalues can be calculated from the matrix equation, 
which follows from (24) 

E J p 2 dp¥ t (p) F {p) = J p 2 dp (cu p + a p ) F t (p) F (p) (34) 

(2n) Jo yUpVLpJo y/oj q Q q 

If the system (24) has been decoupled, or the contribution of nondiagonal elements of the 
IK (p, q) matrix with given radial functions in (34) is zero, Eq. (34) immediately gives the 
perturbative solution for the energy levels. As shown in Ref. [19], the contribution of the 
nondiagonal elements K,^ and iC^i in Eq. (34) to order O (a 4 ) is zero for the ^-mixing states. 
Thus, in the present scheme the energy corrections for ^-mixing states can be calculated 
independently for t = J — 1 and £ = J + 1 states. As a result, all triplet states with 
£ = J =f 1 can be treated as pure states. In the case of spin-mixed states the kernel matrix 
K ope has been diagonalized in the basis of quasi-states (26), however the magnetic part of the 
interaction gives rise to the non-diagonal terms (32). Since we are solving the system (34) 
perturbatively, we can use a new basis \ext) = C[ \sg q )+C' 2 \tr q ) which mixes the quasi-states 
with arbitrary constants C[ and C 2 . This leads to a two-level problem with the solution 

E n ,j,mj = (H u + H 22 )/2± (((H n -H 22 )/2) 2 + H 12 H 21 ) 1/2 , where H n = H% e + Hff, 
H 22 = HZ e + H™\ H 12 = H 21 = H{? = Hlf. In our case \H n - H 22 \ » H 12 H 21 , 
because the difference \H n — H 22 \ is of the order of the fine structure which dominates over 
the hyperfine splitting and the magnetic perturbation H\ 2 . Therefore, we can approximate 

E ni J ym j ~ Hu, H 22 . 

The results are presented in the form 

AE n ^ mj = E n ^ mj - K + ma) + = AE HyJ (« 4 ) + AE™^, (35) 
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where Q 2 = ZQ\. The energy corrections AE n j (a 4 ) due to the kernels K*° pe ) (p,q) were 
obtained previously [19]. The corrections AE riy j (a 4 ) contain spin-spin interactions that lead 
to the HFS which is illustrated in Fig. 1 for the low-lying excited states. A detailed analysis 
of the HFS to O (a 4 ) is provided in [19]. We note that the HFS of the lS 1/2 and 2S 1/2 
states is obtained in agreement with the known Fermi splittings [9], i.e., AE HFS (lS'1/2) = 
(Za) 4 m r (8m r /3M), and AE HFS (25*1/2) = (Za) 4 m r (m r /3M), where M = m\ + m 2 . The 
HFS of states with £ > 0, however, is more complicated [19]. In standard spectroscopic 
notation it has the form 

AE HFS (n, £, s s ) = AE UtJ=e+1 - AE ntJ=ijSa (36) 
(Zafm r 1 (2£+l-£- 1 2m r 1 \ 
" 2 £+ 1 V 4£{£ + l) + M 2£ + 3J ' 



AE HFS (n, £, s t ) = AE ntJ=£tSt - A£„ )J=W (37) 
(Zafm r 1 (2£+ 1 - f" 1 2m r 



n 3 2£+l\ 4£(£ + l) M 2£ — 

where the quantity ^ is defined by Eq. (56), but with the quantum number J replaced by £. 
The formulae (36) and (37) are valid for all quantum numbers n, £ and for any mass values 
mi,m2. The weak external field further splits the energy levels. Eqs. (36) and (37) give 
excellent agreement with experiment for the HFS [19]. 

The energy corrections AE ( f^ lj remove the degeneracy with respect to the mj quantum 
number. The solution of Eq. (34) in the above-made approximation can be written in the 
form of Eq. (1) for all states. 

For all pure states {£ = J =f 1) we obtain the following results: 
for £ = J - 1: 

m 12 J-l / &i, 2 \ 1 fQQ , 
91,2 = 1 j— + I — - - 1 I -, 38 

for 1 = J + l: 

1 mi, 2 J + 2 /g aii3 \ 1 
= 1 - m 1+ m 2 J+l 7TT- (39) 

For spin-mixed states £ = J 7^ the solution of Eq. (34), as mentioned, reduces to a 
standard two-energy level problem. The diagonal elements of the kernel matrix give the 
first-order Zeeman splitting (in O (B)) in the quasi-spin representation (26), which was used 
to derive the HFS energies (36) and (37). Note that the non-diagonal elements give a 
contribution to higher-order Zeeman splitting corrections. 

To first order in the magnetic field strength we obtain the Lande factors to be 



m 1 + m 2 \ 2J(J + 1)J 2 \2J (J + 1) m l +m 2 

mi A _ 1 =F£ \ ( 1 2 \mi -m 2 \ 

m 1 + m 2 \ 2J(J + l)J 2 \2J (J + I) m l +m 2 
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where the upper sign is taken for sg q and lower sign for tr q states respectively. Our ex- 
pressions (40-41) are symmetrical with respect to the masses of the two particles. Ob- 
viously all these first-order Zeeman corrections, AEJ^ , vanish for the positronium case 
(mi = m 2 = m e , Z = 1), as expected. The intrinsic factors g Sl 2 associated with the spins of 
the individual particles can include QED corrections. 

In the case when m 2 >> m\ our general results agree with the result from Eqs. (2,4) 
in which the orbital motion of the heavy particle is ignored. It is only in this limit (as 
discussed below Eq. (27)), that the total angular momenta of the individual particles are 
not related through the common angular momentum £, and can be written as ji — £ ± 1/2, 
and j'2 = 1/2. In J1-J2 coupling, the eigenstates are taken to be the eigenstates of the 

operators jf = + s^j , j| = s 2 ,, J 2 , and J z , and are designated as \ j1j2Jmj) in contrast to 

the spin- mixed \LsJmj) and pure states \LSJmj) which diagonalize the expectation value 
of the Hamiltonian to order 0(a 4 ). To facilitate the comparison we make the following 
replacement of quantum numbers: F — > J, J — > ji, L — > £\ — £, S — > Si, I — > §2- It follows 
that for all pure states £ = J+1, formulae (38-39) and (2-4) give the same result, namely, 

9I = 1+ (^_ 1 )1, 92 ^I, (42) 
for £ = ji - 1/2 = J — 1 and 

9i = 1 -KT- 1 )jTT 1 92 = -TTTT (43) 

for £ = ]i + 1/2 = J + 1. 

In the limit m 2 >> mi the energy levels of spin- mixed states AE e f^ S9q) and AE e f^ trq) 
reduce to AE^ e+1 ^ 2 Jmj and AEj^ e _y 2 Jmj respectively, and the Lande factors given by 
(40-41) take the form 

2J + 3 , f g 8l \ 1 _ 1 ( g S2 



* 2^ If 'J 7- ^ 2 = ~7+T "'2 _1 J 7+T' ^ 

and for £ = ji + 1/2 = J (tr q ) 

2 J - 1 (g Sl \ 1 1 /# S2 \ 1 



£1 = 



2J+ 1 



for £ = ji — 1/2 = J (sg q ). Here the decoupling angle (3 is given by £ « 1/ (2£ + 1) in the 
m 2 >> mi case. 

Formula (4) gives a similar result for the second particle, but for the lighter particle 
Eq. (2) yields 

2^ + 3 (g Sl \ 2J + 3 

5l = 27T7 + lT- 1 J (2j + i)(j + i) (46) 



for (sg q ) states, and 



2J-1 fg Sl ^ 2J-1 



"-wri-KT-ViWTT) (47) 
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for (tr q ) states. This result agrees with (44-45) only in the particular case of g Sl = 2. Note 
that most theoretical and experimental results are concerned with nS\/ 2 (J = 1) states for 
which the "mass ratio" correction in (38) disappears. Thus our results will be most useful 
for I > states. 

In Tables 1-3 we present results of our calculations of the g-factors for the first excited 
states in hydrogen, muonium, and muonic hydrogen respectively. Only states with non- 
zero total angular momentum are included. Eqs. (40-41) are used for the spin-mixed states 
Pi/2(j=i), ^3/2(j=i), £>3/2(j=2), A,/2(j=2), Eq. (38) is used for the pure state P 3 /2(j=2)- 

Table 1. ^-factors for the electron (gi) and proton (g 2 ) respectively in excited atomic 
hydrogen states. Results from the present calculation, Eqs. (38,40) for electrons, are com- 
pared with Eq. (2) in the top half of the table. For protons the bottom half displays the 
present results from Eqs. (38,41) in comparison with Eq. (4). Each row contains in the upper 
part the Lande factor where the intrinsic g s - value is corrected for the anomaly (see text), 
while the numbers below are based upon the Dirac value g s = 2. 



pe 


Pl/2(J=1) 


Ps/2(J=1) 


^3/2(J=2) 


D$/2{J=2) 


£>5/2(J=2) 


g x using Eqs. (38), (40) 


0.33237 
0.33296 


1.66740 
1.66622 


1.00032 
0.99973 


0.59912 
0.59951 


1.40008 
1.39949 


g x using Eq. (2) 


0.33294 
1/3 


1.66765 

5/3 


1.00059 
1 


0.59965 
3/5 


1.39945 
7/5 


g 2 using Eqs. (38), (41) 


1.79321 
1.00036 


-0.89597 
-0.49955 


0.89670 
0.50027 


0.89691 
0.50049 


-0.59711 
-0.33283 


g 2 using Eq. (4) 


1.79285 
1 


-0.89642 
-1/2 


0.89642 
1/2 


0.89642 
1/2 


-0.59762 
-1/3 



Table 2. Same as in Table 1, but for muonium. The Lande factor for the electron is gi, 
and for the muon it is g 2 . 



ji + e 


P\/2(J=l) 


^3/2(J=l) 


^3/2(J=2) 


Dz/2{J=2) 


Db/2{J=2) 


9l using Eqs. (38), (40) 


0.329451 
0.33004 


1.66392 
1.66274 


0.99818 
0.99759 


0.59527 
0.59566 


1.39610 
1.39551 


#i using Eq. (2) 


0.33294 
1/3 


1.66765 

5/3 


1.0006 
1 


0.59965 
3/5 


1.39945 
7/5 


g 2 using Eqs. (38), (41) 


1.00434 
1.00320 


-0.49657 
-0.49598 


0.50299 
0.50241 


0.50491 
0.50433 


-0.32923 
-0.32884 


g 2 using Eq. (4) 


1.00117 
1 


-0.50058 
-1/2 


0.50058 
1/2 


0.50058 
1/2 


-0.33372 
-1/3 



Table 3. Same as in Table 1, but for muonic hydrogen. The Lande factor for the muon 
is gi, and for the proton it is g 2 . 



12 





Pl/2(J=1) 


^3/2(J=l) 


^3/2(J=2) 


^3/2(7=2) 


^5/2(7=2) 


9l using Eqs. (38), (40) 


r\ o *i c\ ^7 

0.26707 


1.58232 


0.95019 


0.50961 


1.30580 


0.26765 


1.58116 


0.94960 


0.51000 


-1 OAT r~\ 1 

1.30521 


g 1 using Eq. (2) 


0.33295 
1/3 


1.66764 

5/3 


1.00058 
1 


0.59965 
3/5 


1.39946 
7/5 


g 2 using Eqs. (38), (41) 


1.85425 


-0.80664 


0.94682 


0.98584 


-0.50225 


1.06317 


-0.41198 


0.55040 


0.58982 


-0.23836 


Q2 using Eq. (4) 


1.79285 


-0.89642 


0.89642 


0.89642 


-0.59762 


1 


-1/2 


1/2 


1/2 


-1/3 



Our calculations (given to five digits after the decimal point) are to be compared with the 
(m 2 — > oo) results (2-4). Upper values for each ^-factor have taken into account the following 
anomalous magnetic moment values: g e /2 = 1.00118, g p /2 = 1.792847 g^/2 = 1.001166 
[1,9,14]. The intrinsic proton anomaly reflects the fact that it is not a fundamental particle, 
while in the case of electrons and muons the lowest-order radiative correction was included. 
The lower values in each row were calculated with g Sl 2 = 2. We used the following values 
for the mass ratios: m p /m e ps 1836.15267 and m^/m e ps 206.76828 [1,9,14]. 

For the case of muonium we find that the deviations between the present results and those 
obtained from the one-body limit are in the few-percent range. The muon as the heavier 
of the two particles acquires a systematically increased Lande factor, while the values are 
always lowered for the electron. 

For muonic hydrogen the effects are more pronounced, and range from 3 to 25 % for the 
states shown in Table 3. Only those results which take the anomalous magnetic moment 
of the proton into account should be considered as physically relevant. The systematics are 
similar to those shown in Table 2 for muonium, with the largest decrease in the Lande factor 
observed for the muon in the P\/2(.j=\) state (-25 %), while the largest increase (19 %) for 
the proton g- value occurs in the D 5 / 2 (j=2) state. 

For atomic hydrogen the effect is smallest due to the small e/p mass ratio. Given that 
atomic spectroscopy is far more advanced in hydrogen than in muonic atoms one should not 
neglect these corrections. For the two above-mentioned states which are most affected we 
observe about 0.1 % deviations in the electron and proton Lande factors respectively. 

As mentioned above, our results are applicable only in low magnetic fields, such that the 
hyperfine energy splitting exceeds the Zeeman splitting, namely 

B « AW "- <,) . (48) 

Thus, formula (48), for 2P 3 / 2 states, requires that B « 300 gauss for muonium and B « 
100 gauss for hydrogen. 

5. Conclusion 

We have used the Hamiltonian variational method in reformulated QED to derive rel- 
ativists stationary-state equations for two-fermion systems in an external magnetic field. 
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These equations can include interactions to any order of the coupling constant, at least 
in principle. The classification of the states follows naturally from the conserved quantum 
numbers which appear in the trial state (6). For given total angular momentum J there are, 
in general, coupled equations, both for mixed-spin states, and for triplet mixed-£ states (cf. 
Eq. (24)). We present explicit forms for the kernels (momentum-space potentials) for the 
case of a constant, weak external magnetic field. 

We solved the radial equations perturbatively to obtain the Zeeman splitting of the HFS 
to order 0(a 4 ), and calculated the (/-factors for the system of two bound fermions. Our 
results are applicable to all states {i.e. for all quantum numbers) and any fermion masses. 
In the limit 1712 » mi our formulae reproduce the well-known </- factor result. For the spin- 
mixed states, however, Eq. (2) is found to be not exact if the intrinsic magnetic moment is 
different from the Dirac value g si =2. 
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Appendix. One-photon exchange kernels for the spin-mixed states 



to order a 4 



We use the notation z = (p 2 + q 2 ) /2pq, and Qj(z) is the Legendre function of the second 
kind [24]. The contributions of the various terms to the kernel are as follows (£ — J (J > 

1), P = (-1) J+1 ): 
(i) orbital term 



2ttQiQ 2 ^ . . 

Mi (P, Q) = 1^22 (P, Q) = Qj( z ) 



+ 



pq 

nQ\Q2 ( f rri! m 2 



(49) 



2mi77l2 V\ m 2 m l 

(ii) spin-orbit interaction 



p q 



— + — -(J-l))[- q + l) Qj(*) + 2 (J + 1) Q J+1 (z) ) , 



/Cir° } (p, q) = 0, 



{P ' q) ~ 2m x m 2 



mi m 2 
m 2 mi 



2^7(7+1) 



^2 °\p <?) ~ 7r< ^ 1< ^ 2 ^ 1711 ' 1712 



(Hi) spin-spin interaction 



.— + — + 4 

2mim 2 \m 2 mi 



K ( i r > (P.9) = 0, 

The diagonalization condition 

tan 2/3 (/C 22 (p, g) - /Cn (p, <?)) = 2/C i2 (p, q) 
determines the parameters (5 and 



2J+ 1 
1 



2J+ 1 



(Qj+i (*) - Qj-i (*)) , 
(Qj+i W-Qj-iW), 



tan 2/5 = 2 



mi — m 2 



mi + m 2 



and 



m l — m 2 

mi + m 2 



-1/2 



J(J+1) + 1 



(50) 
(51) 

(52) 



(53) 



(54) 



(55) 



(56) 
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Therefore, we obtain the diagonalized kernels for the quasi-states 

(57) 



■kQ 1 Q 2 £±1 



2m 1 m 2 £(2J + 1) 



(Qj+i (*) - Qj-i • 
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Fig. 1 . Zeeman splitting of HFS for two-fermion bound state system 



